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Objectives: Lacunes are an important disease feature of cerebral small vessel disease (SVD) but their relation- 
ship to cognitive impairment is not fully understood. To investigate this we determined (1) the relationship 
between lacune count and total lacune volume with cognition, (2) the spatial distribution of lacunes and the 
cognitive impact of lacune location, and (3) the whole brain anatomical covariance associated with these 
strategically located regions of lacune damage. 

Methods: One hundred and twenty one patients with symptomatic lacunar stroke and radiological leukoaraio- 
sis were recruited and multimodal MR1 and neuropsychological data acquired. Lacunes were mapped semi- 
automatically and their volume calculated. Lacune location was automatically determined by projection 
onto atlases, including an atlas which segments the thalamus based on its connectivity to the cortex. Lacune 
locations were correlated with neuropsychological results. Voxel based morphometry was used to create 
anatomical covariance maps for these 'strategic' regions. 

Results: Lacune number and lacune volume were positively associated with worse executive function (num- 
ber p < 0.001; volume p < 0.001) and processing speed (number p < 0.001; volume p < 0.001). Thalamic la- 
cunes, particularly those in regions with connectivity to the prefrontal cortex, were associated with impaired 
processing speed (Bonferroni corrected p = 0.016). Regions of associated anatomical covariance included the 
medial prefrontal, orbitofrontal, anterior insular cortex and the striatum. 

Conclusion: Lacunes are important predictors of cognitive impairment in SVD. We highlight the importance 
of spatial distribution, particularly of anteromedial thalamic lacunes which are associated with impaired 
information processing speed and may mediate cognitive impairment via disruption of connectivity to the 
prefrontal cortex. 

© 2014 Published by Elsevier Inc. 
This is an open access article under the CC BY-NC-ND license 
(http://creativecommons.Org/licenses/by-nc-nd/3.0/). 



1. Introduction 

Cerebral small vessel disease (SVD) is the major cause of vascu- 
lar cognitive impairment in the elderly (Pantoni, 2010) producing a 
characteristic cognitive profile which includes reduced information 
processing speed and executive dysfunction, but relatively preserved 
episodic memory (Charlton et al„ 2006). SVD is thought to impair 
cognition by disrupting the efficiency of white matter pathways that 
connect the networks that underlie cognitive processes (O'Sullivan et 
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al., 2004). Magnetic resonance imaging (MR1) is essential to diagnosis 
and a useful tool to investigate the mechanisms of cognitive impair- 
ment in SVD. A number of MRI markers of cognitive impairment in 
SVD have been proposed, of which T2 white matter hyperintensi- 
ties (WMH) and lacunes (3-15 mm CSF filled cavities of presumed 
vascular origin) are best described (Wardlaw et al., 2013). Here we 
focus on lacunes and the mechanisms by which they lead to cognitive 
impairment in SVD. 

Lacune number has been associated with cognitive impairment 
in some studies in SVD (Gold et al., 2005; Visvanathan et al., 2007; 
Lee et al., 2011), but not in others (Nitkunan et al., 2008). How- 
ever lacunes vary in size and the relationship might be stronger 
with total lacune volume. This has been suggested in CADASIL 
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(Cerebral Autosomal-Dominant Arteriopathy with Subcortical In- 
farcts and Leukoencephalopathy), a young-onset genetic form of SVD 
(Viswanathan et al., 2007), but not yet examined in sporadic SVD. We 
investigate the relationship between both lacune count and total la- 
cune volume with cognition in sporadic SVD and assess whether any 
associations are independent of other pathological markers of SVD 
determined from MR1. 

The impact of lacunes on cognition is likely to depend not only 
on their total load but also on their location. Thalamic lacunes in 
particular have been associated with cognitive impairment in SVD 
(Benisty et al., 2009) and CADASIL (Duering et al., 2011). The thala- 
mus is a complex structure with multiple subcomponents possessing 
unique architecture and connectional properties. Previous studies us- 
ing probabilistic tractography in human brains have showed that the 
thalamus can be segmented into distinct subregions based on cortical 
connectivity patterns that correspond to thalamic nuclei described in 
histological studies (Behrens et al., 2003). In view of this we investi- 
gate the spatial distribution of lacunes to define the specific cognitive 
impact of strategically located lacunes, particularly those in the tha- 
lamus using a thalamic atlas based on its connectivity to the cerebral 
cortex (Johansen-Berg et al., 2005). Guided by the results of our la- 
cune spatial analysis, we identify a region of interest (ROI) within the 
thalamus. We then use voxel based morphometry (VBM) to investi- 
gate whether there are any significant associations between this ROI 
and other brain regions. This technique, termed anatomical (or struc- 
tural) covariance, has previously been used to show that regions that 
are structurally or functionally related 'covary' in grey matter (GM) 
density (Alexander-Bloch et al., 201 3 ; Soriano-Mas et al., 2013). Deter- 
mining 'strategic' lacune locations within the thalamus may improve 
our understanding of the disconnected cortical-subcortical pathways 
underlying cognitive impairment in SVD. 



2. Methods 

2.3. Patients 

Patients with SVD were recruited as part of the prospective 
St George's Cognition and Neuroimaging in Stroke (SCANS) study 
(Lawrence et al., 2013). For this analysis data from the SCANS base- 
line assessment was used. Patients were recruited between March 
2007 and October 2010 from the inpatient and outpatient stroke ser- 
vices of three hospitals covering a geographically contiguous area 
of South London (St George's, King's College and St Thomas' Hospi- 
tals). Inclusion criteria comprised a clinical lacunar stroke syndrome 
(Bamford et al., 1987) with an anatomically corresponding lacunar 
infarct on MRI in addition to confluent white matter hyperintensities 
(WMH) on MRI (Fazekas grade 2 or higher, Fazekas et al., 1993). Ex- 
clusion criteria were as follows: any cause of stroke mechanism other 
than SVD (e.g., cardioembolic source or extra- or intra-cerebral artery 
stenosis of >50%), other major central nervous system disorders, ma- 
jor psychiatric disorders, any other cause of white matter disease, 
contraindications to MRI, or non-fluent in English. The study was 
approved by the local ethics committee and all patients gave writ- 
ten informed consent. MRI acquisitions and cognitive assessments 
were performed at least 3 months after the last stroke to exclude 
acute effects on cognition. Imaging was performed at a median of 
7 months after their symptomatic stroke (minimum: 3 months, max- 
imum: 36 years). 59% of patients were scanned within 12 months 
of their index stroke. All patients were also screened for cardiovas- 
cular risk factors including hypertension (defined as systolic blood 
pressure > 140 mm Hg or diastolic > 90 mm Hg or treatment with 
antihypertensive drugs), hypercholesterolaemia (defined as a serum 
total cholesterol > 5.2 mmol/1 or treatment with a statin), diabetes 
mellitus and smoking. 



2.2. Magnetic resonance imaging acquisition 

Images were acquired on a 1.5 T Signa HDxt scanner (General 
Electric, Milwaukee, WI, USA) with maximum gradient amplitude of 
33 mTm -1 and a proprietary head coil. All image sequences were 
acquired across the whole brain and total imaging time was approx- 
imately 45 min. Patients were placed in the head coil in a neutral 
position with an alignment marker at the nasal bridge to standardise 
head position. Minimal movement was ensured by use of foam pads 
and a Velcro strap across the forehead. 

The imaging protocol included the following: (1) Fluid Attenu- 
ated Inversion Recovery (FLAIR) sequence - TR/TE/TI = 9000/130/ 
2200 ms,field-of-view(FOV) = 240 x 240 mm 2 , matrix =256 x 192, 
28 axial slices of 5 mm thickness and (2) spoiled gradient echo re- 
called Tl -weighted (SPGR) 3D coronal sequence — TR/TE = 11.5/ 
5 ms, FOV = 240 x 240 mm 2 , matrix = 256 x 192, flip angle = 18°, 
176 coronal slices of 1.1 mm thickness reconstructed to an in plane 
resolution of 1.1 mm. 

2.3. Structural image pre-processing 

Structural image preprocessing was performed to compute ac- 
curate deformation field maps to align Tl -weighted volume images 
to a group-average template generated from the data cohort used 
in this study. The preprocessing also involves computation of Tissue 
Probability Maps (TPMs) for grey matter (GM), white matter (WM), 
cerebrospinal fluid (CSF) and white matter hyperintensities (WMH) to 
increase tissue segmentation accuracy across the cohort. Tl -weighted 
images were segmented using New Segment in SPM12 into isotropic 
(1 mm 3 voxel resolution) GM, WM and CSF tissue classes. These were 
used to estimate deformations to a group-average template generated 
from the data cohort using the Shoot toolbox in SPM12 (Ashburner 
and Friston, 2011). The Tl -weighted and FLAIR images were skull- 
stripped using the segmentations at a threshold of 0.1 and warped 
to the group-average template. These were used to create population 
specific TPMs using a modified multivariate mixture of Gaussians 
(mMoG) (Lambert et al., 2013). This was performed to increase seg- 
mentation accuracy across the cohort, and provide TPMs for WMH 
from the FLAIR and Tl weighted images, allowing automatic segmen- 
tation at an individual level. It is necessary to note that GM FLAIR 
hyperintensities were not segmented using this technique. Instead, 
these would have been included in the GM tissue class, and were not 
analysed. Three population specific TPMs for GM, WM and CSF were 
created using the warped Tl -weighted images, and a WMH lesion 
TPM was created from Tl -weighted and FLAIR images. The default 
SPM TPMs were replaced by the population specific TPMs and New 
Segment was re-run on the native space images to generate improved 
segmentation maps for GM, WM and CSF, and a WMH segmentation 
map for each individual. These WMH segmentation maps were bina- 
rised at a threshold set for each individual by checking results man- 
ually to ensure accurate correspondence with lesions on the FLAIR 
image. Results were manually refined where necessary to optimise 
accuracy. 

An additional step was performed to repair the segmentations for 
regions of tissue damage. This was motivated by the observation that 
regions affected by pathology are frequently misclassified during rou- 
tine segmentation. For example, regions of gliosis are misclassified as 
CSF and regions of WMH as GM leading to erroneous deformation 
estimations and subsequent inaccuracies in estimating warps to the 
group-average space. For this reason, WMH regions were automat- 
ically reclassified as WM and any erroneously classified tissue was 
removed from the GM and CSF segmentation maps. Regions with la- 
cunes were also corrected using the most likely tissue type for each 
voxel derived from the population TPMs. All images were visually 
inspected for quality. 

All corrected tissue segmentations (including GM, WM and CSF) 
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were registered to the common group-average isotropic ( 1 mm 3 voxel 
resolution) template using a diffeomorphic registration algorithm 
(Shoot toolbox in SPM 12, Ashburner and Friston, 2011), producing 
a deformation field for each individual to the group-average space. 

2.4. Lacune volume 

Lacunes were identified in native subject space by a consultant 
neuroradiologist, utilising a multimodality view with Tl -weighted, 
T2 -weighted and FLAIR images. A lacune was defined as a CSF filled 
cavity, 3-15 mm in diameter with a surrounding rim of FLAIR hyper- 
intensity (Wardlaw et al., 2013). Cavity size thresholds were chosen 
as lesions that are less than 3 mm in diameter are more likely to be 
perivascular spaces than lacunes and cavities greater than 15 mm 
are less likely to reflect an underlying small vessel disease aetiology 
(Wardlaw et al., 2013). 

To extract lacune regions, the centre voxel of each lacune was iden- 
tified (on Tl -weighted images) and an in-house 6-neighbourhood 
connectivity region-growing algorithm was applied to delineate the 
extent of the lacune. For each subject this algorithm identified a 
threshold boundary for lacune edges (based on the signal intensities 
of all lacune voxels in each brain). Region growing was applied using 
an iterative dilation method and initiated at each central lacune voxel 
until algorithmic termination at the lacune edge. This technique pro- 
vided binary lacune maps for each subject. Lacune maps were visually 
inspected and manually adjusted where this process did not perform 
optimally. Total lacune number (i.e. lacune count) and volume were 
calculated for each subject. 

2.5. Lacune location 

We identified the anatomical location of the already identified 
lacunes with respect to neuroanatomical atlases of (i) white matter 
(WM), (ii) subcortical, and (iii) thalamic structures. To achieve this, 
the previously calculated, population optimised deformation fields 
(see Section 2.3) were used to register the lacune maps to the group- 
average template to create a group-level lacune Tissue Probability 
Map (TPM). Anatomical atlases were used to define the lacune loca- 
tion. These are provided in MNI space, so first needed aligning with 
the group average space. This was done by registering the MNI-152 
Tl -weighted image provided with the FSL-package to the group aver- 
age template using symmetric diffeomorphic non-linear registration 
(Advanced Normalisation Tools, ANTS; http://www.picsl.upenn.edu/ 
ANTS/) (Avants et al., 2011), with parameters optimised for co- 
registration of Tl-weighted images (Klein et al., 2009). The quality of 
the transformation was carefully inspected before these warps were 
applied to three standard atlases provided with the FSL package: 1) 
the label-based JHU DTI white-matter atlas (Mori et al., 2005), 2) 
the probabilistic Harvard-Oxford subcortical structural atlas (http:// 
www.cma.mgh.harvard.edu/), and 3) the Oxford thalamic connectiv- 
ity atlas (Behrens et al., 2003; Johansen-Berg et al., 2005). The Oxford 
thalamic connectivity atlas classifies the thalamus based on its maxi- 
mal connectivity to seven principal cortical regions using probabilis- 
tic tractography, providing "Connectivity Defined Regions" (CDRs). 
To produce non-overlapping, well-defined thalamic CDRs, the atlas 
was thresholded at 25% based on the maximum-probability labels. At 
a threshold of 25%, this atlas has been shown to result in thalamic 
regions corresponding to the nuclei or nuclear groups described in 
previous histological studies (Johansen-Berg et al., 2005). 

For each subject, the lacune volume in each atlas region was com- 
puted. The JHU white matter and thresholded thalamic connectivity 
atlas comprised binary mask images which allowed regional volumes 
to be calculated from overlap with lacune maps in standard space. For 
the probabilistic sub-cortical atlas, the sum of probabilities within 
the lacune voxels was used to calculate the number of equivalent 



whole lacune-voxels, from which the regional lacune volume was 
computed. 

Overlaps between regions and lacunes were visually inspected. 
Lacune overlaps with a volume of less than 10 mm 3 were not reli- 
ably visible within each region. Therefore, for each atlas region, if the 
lacune volume was: (a) less than 10 mm 3 then the region was consid- 
ered to contain no lacunes, but (b) greater than or equal to 10 mm 3 
then the atlas region was considered to contain one or more lacunes. 
This threshold ensured confidence in identified lacune locations and 
ensures that results are representative of a clinically feasible method- 
ology. 

2.6. Brain volume 

Brain volume was estimated from Tl-weighted images us- 
ing the automated SIENAX program (Smith et al., 2002) (FM- 
RIB software library, FSL v4.1; FMRIB analysis group, Oxford, UK; 
http://www.fmrib.ox.ac.uk/fsl ). For each subject, extracted brain tis- 
sue and skull images were registered to a standard brain and skull 
geometry to obtain estimates of both brain volume and brain volume 
normalised with respect to skull size. Errors in tissue classification 
were manually corrected. 

2.7. Lesion load 

White matter hyperintensities on FLAIR were delineated by a sin- 
gle rater using the DISPUNC program (Grimaud et al., 1996; David 
Plummer, University College London, UK). Lesions >2 mm in diame- 
ter were included. Whole brain lesion maps were generated and lesion 
load calculated as a percentage of non-normalised brain volume. 

2.8. Anatomical covariance analysis 

Informed by the findings from the above lacunar spatial analy- 
sis, we performed a further analysis to characterise the anatomical 
covariance associated with the thalamic region of damage. Struc- 
tural images were preprocessed for voxel based morphometry (VBM) 
analysis using Statistical Parametric Mapping software (SPM12, 
2012, http://www.fil.ion.ucl.ac.uk/spm/software/spml2/ ) running 
under MATLAB (201 0b) (Mathworks, Natick, MA, USA). Computed de- 
formation fields (see Section 2.3) were used to warp the corrected CM 
and WM segmentation maps to the group-average template. The re- 
sulting images were modulated by multiplying the correspondingja- 
cobian determinants. Finally the images were smoothed using a 6 mm 
FWHM Gaussian kernel prior to statistical analysis. The preprocessing 
required for this analysis was optimised for our study population as 
described in Section 2.3. 

2.8. 1. Thalamic region of interest 

Informed by the results of our initial lacunar spatial analysis (see 
Results section), only the prefrontal thalamic CDR from the Oxford 
thalamic connectivity atlas was used to create our thalamic region of 
interest. The regions of lacunes that overlapped with the prefrontal 
CDR were selected by multiplying the binarised CDR and lacune TPM 
images. The resulting map was binarised at 1% and used for further 
analysis as the pre-frontal lacunar region of interest (ROI). Prior to 
statistical analysis, the average modulated density of the prefrontal 
lacunar ROI was computed from the smoothed modulated images 
across all subjects. 

2.9. Neuropsychological assessment 

Cognitive assessment was carried out by a neuropsychologist 
within 2 weeks of the MRI scan. Well established standardised tests 
were chosen to include measures sensitive to the pattern of cognitive 
impairment associated with SVD. Premorbid IQwas estimated using 
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the National Adult Reading Test - Restandardised (NART-R) (Nelson 
and Willison, 1991). Tasks were grouped into broad cognitive func- 
tions: executive function, processing speed, working memory, and 
long term (episodic) memory (Lawrence et al., 2013). 

Executive function: Trail making test, (Reitan, 1 955) Modified Wis- 
consin Card Sorting Test (Nagahama et al., 2003), Phonemic Fluency 
(FAS) (Delis et al., 2001). 

Processing speed: Wechsler Adult Intelligence Scale-Ill ( Wechsler, 
1997a) Digit symbol substitution (Wechsler, 1997b), Speed of Infor- 
mation Processing Task (Coughlan et al., 2009), Grooved Pegboard 
Task (Matthews, 1964). 

Working memory: Digit Span Task (Coughlan et al., 2009). 

Long term (episodic) memory: Wechsler Memory Scale-Ill 
(Wechsler, 1997b) Logical Memory and Visual Reproduction. 

Task performance was age scaled using manual published norma- 
tive data, transformed into z-scores and aggregated to construct the 
cognitive indices by averaging across the component test measures 
for each subject. For further details on the cognitive assessment please 
refer to Lawrence et al. (2013). 

The Geriatric Depression Scale (CDS) (Lesher and Berryhill, 1994) 
was used to screen for depressive symptoms and to measure approx- 
imate depression symptom levels. 

2. 1 0. Statistical analysis 

Statistical analysis was performed using Statistical Package for 
Social Sciences version 20.0 (SPSS, Chicago, IL, USA) except for the 
anatomical covariance analysis (Section 2.10.3). Results were con- 
sidered significant at p < 0.05 following Bonferroni adjustment for 
multiple comparisons. All analyses were corrected for age, gender and 
premorbid IQ. Lacune counts and lacune volumes were log 10 trans- 
formed to improve the normality of distributions. 

2.10.1. MRI parameters and cognition 

Partial correlation was used to investigate the relationship be- 
tween cognitive indices and lacune count, lacune volume, WMH vol- 
ume and brain volume. Additional analyses further controlled for hy- 
pertension and for the geriatric depression score as previous stud- 
ies have reported an association between these variables and cogni- 
tive impairment in SVD (Thomas and O'Brien, 2008; Birns and Kalra, 
2009). 

Multiple linear regression analysis was used to investigate 
whether measures of lacune load were associated with cognition. 
In the latter analysis we controlled for WMH volume and normalised 
brain volume. 

2.10.2. Lacune location and cognition 

Multiple linear regression models were used to investigate the 
association between lacune location and each cognitive index. A re- 
gional analysis was conducted for all atlas regions where at least 18 
subjects were affected with a lacune. This threshold was chosen to 
strike a balance between statistical power and region coverage (see 
Supplementary Table SI). Statistical analyses were corrected for ef- 
fects of age, gender, and premorbid IQ. A complete list of investigated 
atlas regions is provided in Supplementary materials S2, S3 and S4. 

2.10.3. Anatomical covariance analysis 

The prefrontal lacunar ROl (see Section 2.8.1 ) was used as the vari- 
able of interest in a multiple regression model in SPM12. Age, gen- 
der and total intracranial volume (TIV) were included as covariates. 
This technique identified the regions of GM or WM that positively or 
negatively correlated with the average GM density of the prefrontal 
lacunar ROL Results were considered significant after Family-Wise 
Error (FWE) correction for multiple comparisons atp < 0.05. 



3. Results 

3.1. Patient demographics 

One patient was excluded from this analysis because of excessive 
movement artefacts in Tl -weighted images. The final sample size was 
120. Demographics are shown in Table 1. 

3.2. Lacune number and distribution 

Cavitated lacunes were present in 99 (83%) of 120 subjects. While 
all patients had clinical lacunar stroke syndrome (with correspond- 
ing MRI lacunar infarction), not all lesions detected on acute diffu- 
sion weighted imaging subsequently cavitate on Tl -weighted images 
(Loos et al., 2012). The anatomical location of lacunes with reference 
to the subcortical and white matter standard space atlases is shown 
in Fig. 1. Lacunes were most frequently identified within the thala- 
mus (see Table 4). Lacune count was positively associated with WMH 
volume. The association remained significant after controlling for age 
and cardiovascular risk factors (r = 0.423, p < 0.001 ). 

3.3. Association of MRI parameters with cognition (Table 2 ) 

Associations between MRI parameters and cognitive scores con- 
trolling for age, gender and premorbid IQ were strongest with execu- 
tive function and processing speed. Lacune count and lacune volume 
showed negative associations of similar magnitude, although for la- 
cune count the partial correlation coefficients were slightly greater. 
For lacune count and volume there were weaker associations with 
working and episodic memory. Brain volume was strongly associated 
with all cognitive domains. In contrast, associations between WMH 
volume and cognition were weaker. Additional analyses controlling 
for hypertension and depression did not affect the significance of our 
results. 

The relationships of lacune count and lacune volume with exec- 
utive function and processing speed remained significant after con- 
trolling for brain volume and WMH volume (Table 3). 

3.4. Associations between lacune location and cognition 

Regional analysis was performed on the following regions: subcor- 
tical grey matter regions which included the caudate, thalamus and 
putamen and white matter regions which included the internal cap- 
sule, external capsule, superior longitudinal fasciculus, and anterior, 
superior and posterior corona radiata. Thalamic CDRs included the 
prefrontal, premotor, temporal, primary motor and posterior parietal 
cortex CDRs (Fig. 2). 

In subcortical GM, thalamic lacunes were associated with impaired 
processing speed (Table 4; p < 0.001 ), in WM regions, there were asso- 
ciations between lacunes in the superior corona radiata and impaired 
processing speed (p < 0.035), and between lacunes in the posterior 
internal capsule and impaired executive function (p < 0.001 ). The as- 
sociation between thalamic lacunes and impaired processing speed 
survived Bonferroni correction (n = 9, corrected p = 0.027; Table 4). 

The relationship between thalamic connections to specific cortical 
regions and impaired processing speed was further explored (Fig. 2 
and Table 4). Associations were present between impaired processing 
speed and lacunes in the thalamic CDRs with connections to the pre- 
frontal cortex (p < 0.004), posterior parietal cortex (p < 0.016), and 
temporal cortex (p < 0.01 1 ). The association for lacunes in the thala- 
mic CDR with connections to the prefrontal cortex survived Bonfer- 
roni correction for multiple comparisons (n = 4, corrected p = 0.016). 
In order to confirm this we performed a supplementary step by step 
ascending multiple regression analysis where only the association 
with the prefrontal CDR and processing speed remained significant. 
This localised finding motivated the following VBM analysis. 
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Table 1 

Patient demographics and mean MRI parameters. Hypertension was defined as systolic blood pressure > 140 mm Hg or diastolic > 90 mm Hg or those on antihypertensive 
treatment. Hypercholesterolaemia was defined as a serum total cholesterol > 5.2 mmol/1 or treatment with a statin. 



Demographics and risk factors 




N =120 


Mean age (SD) 




70(9.8) 


Mini mental test score (mean (SD)) 




27.6 (2.7) 


Gender 


Female 


42 (35.0%) 




Male 


78 (65.0%) 


Hypertension 


No 


9 (7.5%) 




Yes 


111 (92.5%) 


Hypercholesterolaemia 


No 


17(14.2%) 




Yes 


103(85.8%) 


Smoker 


Never 


55(45.8%) 




Current 


23 (19.2%) 




Ex-smoker 


42 (35.0%) 


Treated diabetes mellitus 


Yes 


22(18.3%) 




No 


96 (80.0%) 




Untreated 


2(1.6%) 


Rankin disability score 


0 


38 




1 


48 




2 


15 




3 


16 




4 


3 


MRI measures at baseline (per subject) 






Mean (SD) number of lacunes 




4.18(5.44) 


Mean (SD) lacune load as % of total brain volume 




0.0754(0.0979) 


Mean (SD) WMH load as % of total normalised brain volume 


3.16(2.63) 


Mean (SD) normalised brain volume in ml 




1295.64(91.29) 



Table 2 

Partial correlation between MRI parameters and cognition: R (p values) controlling for age, gender, and NART IQ. 



Episodic 

Executive function Processing speed Working memory memory 



Lacune count -0.377(0.001) -0.430(0.001) -0.207(0.028) -0.189(0.045) 

Lacunevolume -0.310(0.001) -0.388(0.001) -0.236(0.012) -0.255(0.006) 

WMH volume -0.104(0.274) -0.239(0.011) -0.082(0.389) -0.135(0.153) 

Brain volume 0.328(0.000) 0.403(0.000) 0.285(0.002) 0.281 (0.003) 



3.5. Anatomical covariance: grey matter 

Several grey matter regions positively correlated with the average 
density of the prefrontal lacunar ROl (Fig. 3) which was within the 
prefrontal thalamic CDR. Subcortically, these included the entire tha- 
lamus, caudate, putamen and globus pallidus. Bilateral cortical asso- 
ciations were found in the orbitofrontal cortex, the medial prefrontal 
cortex, the dorsoanterior insular cortex, the parieto-occipital junction 
and the cuneus. No significant negative correlates were found. 

3.6. Anatomical covariance: white matter 

Predominately anterior white matter regions positively co-varied 
with the average density of the prefrontal lacunar ROI (Fig. 4). These 
included the anterior thalamic radiation, the anterior and superior 
corona radiata, the genu and body of the corpus callosum and the 
internal capsule. No significant negative correlates were found. 

4. Discussion 

Our results show that both the number and the volume of lacunes 
are independent predictors of impaired processing speed and exec- 
utive function in symptomatic small vessel disease. We have shown 
that lacunes in the thalamus are independently associated with the 
processing speed deficits. These were further characterised using Ox- 
ford's thalamic connectivity atlas (Johansen-Berg et al., 2005), iden- 
tifying the region with high connectivity to the prefrontal cortex, 
located in the anteromedial thalamus, as significant. We provided 
support for this finding by examining the anatomical covariance as- 
sociated with the lacunar ROl within this region, demonstrating sig- 
nificant regions predominately in the prefrontal, orbitofrontal and 



anterior insula cortex. 

Previous studies have suggested an association between lacunes 
and impaired cognition (Vermeer et al., 2007), but none so far 
have shown this in sporadic symptomatic small vessel disease. The 
Leukoaraiosis and Disability (LAD1S) study, which included patients 
with leukoaraiosis from different clinical settings, reported that la- 
cune count predicted worse executive function (Benisty et al., 2009). 
These results are in agreement with postmortem data (Gold et al., 
2005). In CADASIL the number of lacunes was found to be the most 
significant MRI parameter associated with cognitive impairment (Lee 
et al., 2011). However not all studies have shown such associations, 
including one study in patients with symptomatic lacunar infarction 
(Nitkunan et al., 2008). Our data suggest that lacunes are an impor- 
tant predictor of cognitive impairment in sporadic symptomatic SVD, 
in a similar fashion to their role in the genetic form of the disease 
CADASIL. One previous study in CADASIL found that total lacune vol- 
ume was more strongly associated with cognition than lacune number 
(Visvanathan et al., 2007). However, we found no significant differ- 
ence between the two parameters in sporadic SVD; indeed associ- 
ations with executive function and processing speed tended to be 
stronger for lacune count. In contrast we found weak or absent as- 
sociations between WMH lesion volume and cognition. This may be 
because high signal on T2 -weighted imaging fails to distinguish be- 
tween minor and severe white matter damage. It is also possible that 
by selecting patients with confluent WMH in our inclusion criteria we 
have limited the range of WMH volume values and this could have 
reduced the power to detect a relationship between WMH and cogni- 
tion. However, despite this inclusion criterion there was a wide range 
of WMH volumes from 0.145 to 12.6% of total brain volume. 

Differences between the results of previous studies correlating 
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Fig. 1. Lacune lesion location map showing spatial distribution throughout the brain. The map is superimposed onto the group-average Tl -weighted 1 mm isotropic template and 
shown using the neurological viewing convention. Row A — atlas brain regions illustrating subcortical structures and white matter regions (standard atlases provided by the FSL 
package). A colour key to atlas regions is shown at the bottom of the figure. Row B — distribution of lacunes in MNI space. The colour bar at the bottom right indicates the number 
of lacunes present at each voxel. Row C — lacunes shown superimposed on atlas locations. 

Table 3 

Associations between lacune count, volume and cognition (multiple linear regression analysis controlling for age. gender. NART, brain volume and WMH volume). 

MRI parameters Executive function Processing speed 

fi Significance |3 Significance 

Lacune count -0.222 0.005 -0.261 0.002 

Lacune volume -0.172 0.031 -0.228 0.006 



MRI parameters with cognition may also reflect the concept that SVD 
may have multiple pathologies. It has been suggested that there are at 
least two distinct pathologies underlying SVD; microatheroma caus- 
ing larger, usually isolated lacunar infarcts and a diffuse arteriopathy 
characterised by lipohyalinosis associated with multiple smaller lacu- 
nar infarcts and leukoaraiosis (Fisher, 1968; Boiten et al., 1996). This 
distinction is supported by pathological, radiological and risk factor 
data (Khan et al., 2007). In view of these possible differences and to 
study a homogeneous group of patients we recruited only patients 
with symptomatic lacunar infarction confirmed on MRI as well as 
confluent WMH. We feel this allows a more robust determination of 
the effect of lacunes in a population where the cause of both WMH and 
cognitive impairment is likely to be ischaemic. We recognise that this 
may limit the generalisability of the findings. However the findings 



remain relevant to a large number of stroke patients — about 20-25% 
of all ischaemic stroke is lacunar due to SVD and of these cases about 
half fall into the category of lacunar stroke and WMH. 

The cognitive impairment in SVD associated with lacunes may 
well be due to damage to cortical-subcortical pathways, disrupt- 
ing the complex and distributed networks that underpin processes 
such as executive function and information processing (Bressler and 
Menon, 2010). If this were the case then the spatial location of dam- 
age would be more important than lesion load in determining cog- 
nitive outcome. This work demonstrates that thalamic lacunes are 
associated with cognitive impairment, consistent with previous ob- 
servations (Benisty et al., 2009; Duering et al., 2011). However, the 
thalamus is a complex structure consisting of multiple subnuclei that 
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Fig. 2. Oxford thalamic connectivity probability atlas superimposed on to group-average Tl-weighted 1 mm isotropic template and shown using the neurological viewing 
convention. The colour key at the bottom of the figure represents the classification of thalamic voxels according to the cortical region with which they are connected with highest 
probability. Row A — thalamic regions according to the cortical region with which they are connected with highest probability. The colour key at the right of the figure represents 
this classification. Row B — lacunes shown superimposed on thalamic atlas locations. The colour bar on the right indicates number of lacunes present at each voxel. Row C — axial 
thalamic section from a cytoarchitectonic atlas (Morel et al., 1997). The nuclei have been coloured according to their major cortical connection site. 




Fig. 3. VBM analysis showing anatomical covariance maps of grey matter regions that positively co-varied with the prefrontal lacunar ROl average grey matter density. Statistical 
maps are shown superimposed on to the common group-average 1 mm isotropic template. The colour bar represents the t-statistic. Images are displayed in axial (row 1 ), sagittal 
(row 2) and coronal (row 3) orientations using the neurological viewing convention. 
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Fig. 4. VBM analysis showing anatomical covariance maps of white matter regions that positively co-varied with prefrontal lacunar ROI average grey matter density. Statistical 
maps are shown superimposed on to the common group-average 1 mm isotropic template. The colour bar represents the t-statistic. Images are displayed in axial (row 1 ), sagittal 
(row 2) and coronal (row 3) orientations using the neurological viewing convention. 



Table 4 

Spatial distribution of lacunes and associations with cognition (linear regression analysis controlling for age. gender and NART). (* represents regions that survived Bonferroni 
correction). 

Number of subjects with 
lacunes (% of total 



Region 


sample) 


Processing speed 




Executive function 








P 


Significance 


P 


Significance 


Subcortical grey matter 












Thalamus (whole) 


47 (39.2) 


-0.259 


0.001* 


-0.142 


0.047 


Putamen 


33 (27.5) 


-0.116 


0.153 


-0.104 


0.151 


Caudate 


35 (29.2) 


-0.080 


0.314 


-0.140 


0.049 


White matter structures 












Anterior limb of 


28 (23.3) 


-0.107 


0.174 


-0.046 


0.520 


internal capsule 












Posterior limb of 


28 (23.3) 


-0.146 


0.065 


-0.242 


0.001 * 


internal capsule 












Anterior corona radiata 


18(15.0) 


-0.148 


0.073 


-0.068 


0.362 


Superior corona 


39 (32.5) 


-0.168 


0.035 


-0.110 


0.124 


radiata 












Posterior corona 


19(15.8) 


-0.115 


0.145 


-0.181 


0.010 


radiata 












Superior longitudinal 


21 (17.5) 


-0.024 


0.771 


-0.083 


0.258 


fasciculus 












External capsule 


32 (26.7) 


-0.158 


0.053 


-0.143 


0.052 


Thalamic connectivity 












Prefrontal 


44(36.7) 


-0.225 


0.004 * 


-0.186 


0.046 


Premotor 


23(19.1) 


-0.132 


0.096 


-0.075 


0.291 


Posterior parietal 


37(30.8) 


-0.190 


0.016 


-0.108 


0.127 


Temporal 


28 (23.3) 


-0.204 


0.011* 


-0.174 


0.016 
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have unique functional properties. We demonstrate a specific asso- 
ciation between the anteromedial thalamus and impaired process- 
ing speed. To further validate this in our cohort, we characterised 
the anatomical covariance associated with the area damaged within 
this region, which corresponded well with the described connectiv- 
ity. Contained within this region are the ventral anterior (VAN) and 
the mediodorsal (MDN) nuclei, (Johansen-Berg et al., 2005) (Fig. 2C). 
Tract tracing studies have shown that this region of the thalamus 
has extensive connections to the orbitofrontal, prefrontal and cingu- 
late cortices (Ray and Price, 1992; Klein et al., 2010), which has been 
replicated in-vivo using probabilistic diffusion tractography (Behrens 
et al, 2003; Johansen-Berg et al., 2005; Klein et al., 2010). In addition 
cytoarchitectonic studies and functional neuroimaging studies have 
also described connections between the mediodorsal and ventral nu- 
clei of the thalamus and the anterior insular cortex (Cho et al., 2012; 
Metzger et al., 2013). The combined evidence therefore suggests that 
lacunes in this region may impair processing speed by disrupting 
connectivity to the orbitofrontal and prefrontal cortices. 

The thalamus is also densely interconnected with the striatum and 
cortex via parallel corticostriatal-thalamo-cortical loops (Draganski 
et al., 2008; Metzger et al., 2013). Initially described by Alexander et 
al. (1986), these are arranged topographically into distinct functional 
circuits. The thalamic region identified in this work would overlap 
with the previously identified orbitofrontal and medial-prefrontal 
networks. The structural covariance from our work corresponds well 
with the described connectivity of these functional circuits (Draganski 
et al., 2008), providing indirect supportive evidence for their involve- 
ment. Previous studies in multiple sclerosis (Batista et al., 2012) and 
ageing (Van Der Werf et al., 2001) have reported that thalamic atro- 
phy is associated with impaired processing speed. Our results provide 
a sound biological basis for these observations. It is likely that pro- 
cessing speed is not due to a single neural system, but rather is a 
reflection of activity across several overlapping networks. It may be 
that the anteromedial thalamus forms a 'hub' for the networks in- 
volved in processing speed performance, with impairment resulting 
from a lacune mediated disconnection syndrome. 

Some limitations of our study are that some functionally relevant 
locations of lacunes may have been missed due to limited overlap of 
lacune locations between subjects and due to the conservative nature 
of Bonferroni correction. Furthermore, although the atlas registration 
was manually checked for quality, a perfect image co-registration can- 
not be guaranteed by application of current techniques. The analysis 
of grey matter lacune location is likely to be more robust than that 
for white matter regions which are less well defined than grey mat- 
ter structures. The thalamus had the highest frequency of lacunes and 
therefore associations in this region may have been easier to detect. It 
is possible that lacunes in other locations are also strategically impor- 
tant but we did not have a sufficient sample size to detect significant 
associations. In addition we lacked sufficient numbers to investigate 
the impact of multiple lacunes in one region and interactions between 
lesion locations. A further limitation is that the biological basis for 
structural covariance is incompletely understood and does not nec- 
essarily infer structural connectivity. Structural covariance can also 
be influenced by developmental, genetic and environmental factors. 
Therefore we supported our observations with previously reported 
white matter connectivity results from both healthy in vivo popula- 
tions and ex vivo histology. Despite its limitations, because structural 
covariance does not depend upon intact white matter pathways, it 
may allow the exploration of network hypotheses, particularly in dis- 
eased brains such as in patients with SVD, as it avoids the method- 
ological issues that arise from attempting to use in vivo tractography 
to define networks of interest across regions of degenerating white 
matter tracts. 



5. Conclusion 

Lacunes are an important marker of impaired information process- 
ing speed and executive function, which are characteristic of vascular 
cognitive impairment in SVD. We demonstrated that thalamic lacunes 
impair processing speed, and informed by this observation further 
localised it to the anteromedial thalamus. Lacunes in this area may 
impair processing speed by disrupting connectivity to the prefrontal 
cortex. This study highlights the importance of the spatial distribu- 
tion of lacunes and provides support for the role of corticostriatal- 
thalamo-cortical loops in processing speed deficits in SVD. 
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